Level Crossing Analysis of Burgers Equation in 1 + 1 Dimensions 
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We investigate the average frequency of positive slope , crossing the velocity field u(x) — u = a in 
the Burgers equation. The level crossing analysis in the inviscid limit and total number of positive 
crossing of velocity field before creation of singularities are given. The main goal of this paper is to 
show that this quantity, v& , is a good measure for the fluctuations of velocity fields in the Burgers 
turbulence. 
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I. INTRODUCTION 

The Burgers equation is the simplest nonlinear gener- 
alization of the diffusion equation. The N-dimensional 
forced Burgers equation 



d t u + (u. V)u = vV 2 u + f(x, t) 



(1) 



which describes the dynamics of a stirred, pressure less 
and vorticity-free fluid, has found interesting applications 
in a wide range of non-equilibrium statistical physics 
problems. It arises, for instance, in cosmology where it 
is known as the adhesion model [1], in vehicular traffic 
[2] or in the study of directed polymers in random media 
[3] . In the Burgers equation if the velocity field is a gra- 
dient field u(x, t) = —Vip(x, t) and the random force is a 
gradient random force f(x,t) = —VF(x,t), then associ- 
ated Hamilton Jacobi equation, satisfies in the following 
equation as: 



W = ^ + J(W) 2 + F(f, t) 



(2) 



where v is the viscosity, recently it has been frequently 
studied as a nonlinear model for the motion of an in- 
terface under deposition [4]. The case with large-scale 
forcing was considered in Refs. [5,6] as a natural way 
to pump energy in order to maintain a statistical steady 
state. Burgers equation is then a simple model for study- 
ing the influence of well- understood structures (shocks, 
preshocks, etc) on the statistical properties of the flow. 
As it is well known, eq.(l) in the limit of vanishing viscos- 
ity (y — > 0) displays after a finite time dissipative singu- 
larities, namely shocks, corresponding to discontinuities 
in the velocity field. In the presence of large-scale forc- 
ing, it was recently stressed for the one-dimensional case 
[7,8] and also for higher dimensions [9-11], that the global 
topological structure of such singularities is strongly re- 
lated to the boundary conditions associated to the equa- 
tion. More precisely, when, for instance, space period- 
icity is assumed, a generic topological shock structure 



can be outlined. It plays an essential role in understand- 
ing the qualitative features of the statistically station- 
ary regime. So far, the singular structure of the forced 
Burgers equation was mostly investigated in the case of 
finite-size systems with periodic boundary conditions. It 
is however frequently of physical interest to investigate 
instances where the size of the domain is much larger 
than the scale, so as to examine, for example, the role of 
Galilean invariance [12]. 

Here we describe the level crossing analysis in the context 
of vorticity - free fluid. In the level crossing analysis we 
are interested in determining the average frequency ( in 
spatial dimension ) of observing of the definite value for 
velocity function u(x) — u = a in fluid, v+ , from which 
one can find the averaged number of crossing the given 
velocity in sample with size L. The average number of vis- 
iting the velocity u(x) — u — a with positive slope will be 
= v^L. It can be shown that the i/+ can be written 
in terms of joint probability distribution function (PDF) 
of u{x) — u and its gradient. Therefore the quantity v+ 
carry the whole information of fluid which lies in joint 
PDF of velocity and its gradient fluctuations. This work 
aims to study the frequency of positive slope crossing (i.e. 
v£) in time t on the vorticity - free fluid in a sample with 
size L. We describe a quantity N^ ot which is defined as 
N+ ? +co - 

Jv tot 



f °° v+da to measure the total number of cross- 

J — oo a 

The N+ t 



ing the velocity of fluid with positive slope 
and the path which is constructed velocity of fluid are 
in the same order. It is expected that in the station- 
ary state the N^~ ot to become size dependent. Although 
we exactly determine the velocity dependence of v£ for 
Burgers equation in the inviscid limit and before creation 
of singularities, we compute the time dependence of N^ ot 
(z/+) numerically 

This paper is organized as follows: In section II we dis- 
cuss the connection between z/+ and underlying prob- 
ability distribution functions (PDF) of a fluid [13]. In 
section III we derive the integral representation of z/+ for 
the Burgers equation in 1+1 dimensions and in the invis- 
cid limit before the creation of singularities. Section IV 



closes with a discussion of the present results. 



II. THE LEVEL CROSSING ANALYSIS OF 
STOCHASTIC PROCESSES 

Consider a sample function of an ensemble of functions 
which make up the homogeneous random process u(x,t). 
Let ra+ denote the number of positive slope crossing of 
u(x) — u = a in time t for a typical sample size L (see 
fig.(l) ) and let the mean value for all the samples be 
N+(L) where: 



N+(L)=E[n+(L)]. 



(3) 



Since the process is homogeneous, if we take a second 
interval of L immediately following the first we shall ob- 
tain the same result, and for the two intervals together 
we shall therefore obtain: 



N+(2L) = 2N+(L), 



(4) 



from which it follows that, for a homogeneous process, 
the average number of crossing is proportional to the 
space interval L. Hence: 



or: 



N+(L) cx L, 



N+(L) = u+L. 



(5) 



(6) 



which is the average frequency of positive slope cross- 
ing of the level u{x) — u = a. We now consider how the 
frequency parameter i/+ can be deduced from the under- 
lying probability distributions for u{x) — u. Consider a 
small length dl of a typical sample function. Since we 
are assuming that the process u(x) — u is a smooth func- 
tion of x, with no sudden ups and downs, if dl is small 
enough, the sample can only cross u(x) — u = a with 
positive slope if u{x) — u < a at the beginning of the in- 
terval location x. Furthermore there is a minimum slope 
at position x if the level u(x) — u — a is to be crossed in 
interval dl depending on the value of u(x) — u at location 
x. So there will be a positive crossing of u(x) — u = a in 
the next space interval dl if, at position x, 



u(x) — u < a and 



d(u(x) 



dl 



> 



a - 



(u(x) 



dl 



(7) 



Actually what we really mean is that there will be high 
probability of a crossing in interval dl if these conditions 
are satisfied [14,15]. 

In order to determine whether the above conditions are 
satisfied at any arbitrary location x, we must find how 
the values of y — u(x) — u and y' = ^ are distributed by 




■ h = 



FIG. 1. positive slope crossing of the level u(x) — u = a. 



considering their joint probability density p(y,y'). Sup- 
pose that the level y = a and interval dl are specified. 
Then we are only interested in values of y < a and values 
of y' = (^) > ^jf-, which means that the region between 
the lines y = a and y' = in the plane (y, y'). Hence 
the probability of positive slope crossing of y = a in dl 
is: 



f 

Jo 



dy 



P CX. 

' I dyp(y,y'). 

J a — y'dl 



When <iZ — *• 0, it is legitimate to put: 

p{y,y') =p(y = a,y')- 



(8) 



(9) 



Since at large values of y and y' the probability density 
function approaches zero fast enough, therefore cq.(8) 
may be written as: 



dy' 



dyp(y = a,y'), 



(10) 



a — y'dl 



in which the integrand is no longer a function of y so 
that the first integral is just: J^_ y , dl dyp(y — a,y') = 
p(y = a,y')y'dl, so that the probability of slope crossing 
of y — a in dl is equal to: 



dl / p(a,y')y'dy', 
Jo 



(11) 



in which the term p(a, y') is the joint probability density 
p(y, y') evaluated at y = a. 

We have said that the average number of positive slope 
crossing in scale L is v+L, according to eq.(6). The av- 
erage number of crossing in interval dl if therefore v+dl. 
So average number of positive crossing of y = a in inter- 
val dl is equal to the probability of positive crossing of 
y = a in dl, which is only true because dl is small and 
the process y(x) is smooth so that there cannot be more 
than one crossing of y = a in space interval dl, Therefore 
we have v^dl = dl p(a,y')y'dy' , from which we get 
the following result for the frequency parameter in 
terms of the joint probability density function p(y, y') as 
follows: 



F 

Jo 



p(a,y')y'dy'. 



(12) 



In the following section we are going to derive the v+ 
via the joint PDF of u(x) — u and velocity gradient. To 
derive the joint PDF we use the master equation method 
[16]. This method enables us to find the in terms of 
generating function. 



III. FREQUENCY OF A DEFINITE VELOCITY 
WITH POSITIVE SLOPE FOR BURGERS 
EQUATION BEFORE THE SINGULARITY 
FORMATION 

As mentioned in section I, in the presence of a random 
force f(x, t) , the velocity field of burgers equation in 1 + 1 
dimensions evolves as: 



du du d 2 u ,. , 



(13) 



Differentiating the Burgers equation (eq.13) respect to x, 
we have: 



duj duj 



d 2 



^ + u S + " 2 = i^ + / * (M) ' 



(14) 



where v > 0, u — |^ and f(x, t) is a random force, with 
a Gaussian distribution of mean zero and second moment 
given by: 

(fix, t)f(x', t')) = 2D D(x - x')S(t - t') (15) 

where D(x) is space correlation function and is an even 
function of its argument. It has the following form: 



D(x-x') = 



(16) 



where a is the standard deviation of D(x — x'). Typically 
in the realistic problem, the correlation of forcing is con- 
sidered as smooth function for mimicking the long range 
correlation. We regularize the smooth function correla- 
tion by a gaussian function. When the variance a is in 
order of the system size, we would expect that the model 
would represent a long range character for the forcing. 
So we should stress that our calculations are done for fi- 
nite a ~ L, where L is the system size. Parameters v 
and Do are describing kinematics viscosity and the noise 
strength, respectively. Once trying to develop the sta- 
tistical theory of the Burgers equation it becomes clear 
that the inter-dependency of the velocity field and veloc- 
ity gradient statistics would be taken into account. The 
very existence of the non-linear term in the Burgers equa- 
tion leads to development of the singularities in a finite 
time and in the inviscid limit i.e. v — > 0. So one would 
distinguish between different time regimes. Recently it 
has been shown that starting from the flat interface the 



Burgers equation will develop shock singularity after time 
scale t* , where t* depends on the forcing properties as 
t* ~ AT 1/3 cr [11,16,17]. This means that for time scales 
before t* the relaxation contribution tends to zero when 
v — > 0. In this regime one can observe that the generat- 
ing function equation is closed [11,17-19]. Let us define 
the generating function Z(X, fi, x, t) as: 



Z(X,n,t) = (e 



-i\(u(x ,t) —u) —ifj,uj(x ,t) 



> = (©>• (17) 



Assuming statistical homogeneity i.e. Z x = it follows 
from equations (13) and (14) that Z satisfies in the fol- 
lowing equation: 



dz_ 

~dt 



(18) 



Using the Novikove theorem [11] gives: 
dZ 



dt 



= iZ^ - H iZ^ - X A k(0)Z + n A k xx (0)Z (19) 



where k(x - x 1 ) = 2D D(x - x'), fc(0) = and 
k X x(0) = — ^§° 3 ■ The solution of eq.(19) by using the 
separation of variables is as follows: 

Z(t,n,\) -- 



k(0)t 



e -y m t x ZiM 



(20) 



where Zi(fi,t) satisfies in the following equation: 
dZ 

-j£= iZ llM - iliZ x ^ - \ 2 k{Q)Z 1 
+ [^ 2 k xx (0) + C(t)}Z 1 



(21) 



C(t) is an arbitrary function which should be determined 
by initial conditions. 

The joint probability density function of u and id can 
be obtain by Fourier transform of the generating func- 
tion: 

P{u,uj,t) = -!- / d\d^ u -^ +l ^Z(X,^t), (22) 

27T J 

so by Fourier transforming of the eq.(19) we get the 
Fokkcr-Planck equation as: 

ft ftP ft 2 P ft 2 P 

-P = 3w P + w »_ + fc (0)^-M0)^. (23) 

The solution of the above equation can be separated as 
P(u,w,t) = pi(u, t)p2(cj, t) (for motivation see [20]). Us- 
ing the initial conditions P\(u,0) = S(u) and P2(a>,0) = 
S(uj) it can be shown that: 



P(u, u>, t) 



=e 4fc <o)t x p 2 {u>,t) 



y/4irk(0)t 

where p 2 (LJ,t) is a solution of the following equation: 



(24) 




FIG. 2. plot of u£ vs a for the Burgers equation in the 
strong coupling and before the creation of singularity for time 
scale t/t* = 0.02,0.03 and 0.04. 
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FIG. 3. log-log plot of Nt ot vs t for the Burgers equation in 
the strong coupling before the creation of singularity in the 
velocity field. 



ik p2 = w iB ~ Kx (o) + [3u} + G{t)]p2 (w ' t] 

(25) 

G(t) is an arbitrary function which should be determined 
by initial conditions. Up to now we obtained that, for ev- 
ery time scale, the level cross v£ has a Gaussian behavior 
in terms of a, now for determination of time dependence 
of i/+ we have to determine p2(u>,t). Since the eq.(25) 
is so complex, we solve it using the numerical methods 
[21]. The frequency of repeating a definite velocity field 
(u(x) — u — a) with positive slope can be calculated as 
v a = Jo ojP{a,uj,t)dui. Fig. 2 shows z/+ for various time 
scales before creation of singularity. To derive the N^~ ot 
let us express N^ ot as: 

/ + OG />OC 
da I cuP(a,uj,t)duj (26) 
-oo JO 

Using the numerical integration of eq.(26) one finds 
N+ t ~ t where /3 = 0.50 ± 0.01, In fig. (3) we plot the 
N^ ot as a function of t. 

IV. CONCLUSION 

We obtained some results in the problems of Burg- 
ers equation in 1+1 dimensions with a Gaussian forcing 
which is white in time and Gaussian correlated in space, 
typically in the physical case. We determined the average 
frequency of crossing i.e. v+ of observing of the definite 



value for velocity field u — u — a, from which one can find 
the averaged number of crossing the given velocity field 
in a sample with size L. The integral representation of v£ 
was given for the Burgers equation in the inviscid limit 
before the creation of singularity and it was shown that 
the velocity dependence of the v+ is Gaussian. We apply 
the quantity N^ ot = v+da, which measures the to- 
tal number of positive crossing of velocity and show that 
for the Burgers equation in the inviscid limit and before 
the creation of singularities N* ot scales as t 1 / 2 . 
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